Magnetic-induced phonon anisotropy in ZnCr 2 04 from first principles 
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We have studied the influence of magnetic order on the optical phonons of the geometrically 
frustrated spinel ZnCr2 04 from first-principles. By mapping the first-principles phonon calculations 
onto a Heisenberg-like model, we developed a method to calculate exchange derivatives and subse- 
quently the spin-phonon couping parameter from first-principles. All calculations were performed 
within LSDA+U. 
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The interplay between spin and lattice degrees of free- 
dom give rise to an amazing variety of different phe- 
nomena, including the spin- Teller and the magneto- 
electric |2J effects. Perhaps the simplest example of a 
spin-lattice coupling is the influence of magnetic order 
on phonons. Though first-principles density functional 
methods have been highly successful in describing the 
structural and magnetic properties of dielectrics [3| , lit- 
tle has been done to address the coupling of spins and 
phonons. fn this letter we present for the first time 
an approach to predict the influence of magnetic order 
on optical phonons from first-principles and apply this 
approach to the strongly geometrically frustrated spinel 
zinc chromite. Our method is general and has also been 
recently applied to provide a new understanding of the 
phonon anomalies observed at the ferromagnetic (FM) 
transition in the chalcogenide spinels p|. 

Zinc chromite, ZnCr204, crystallizes in a spinel struc- 
ture. The Cr 3+ ions (S= 3/2) form a network of vertex- 
sharing tetrahedra (a pyrochlore lattice) with strong anti- 
ferromagnetic (AFM) interactions between nearest neigh- 
bor spins, 6cw~ 400K. As a result, zinc chromite is 
strongly geometrically frustrated as evidenced by the 
rather low T^ = 12. 5K (compared to Qcw)- The 
physics of ZnCi'204 at Tjy involves a first-order cubic-to- 
tetragonal (c-t) structural transition (to relieve geomet- 
ric frustration) as it enters the Neel state, bypassing a 
bond-ordered state Q. Recently, Sushkov et al. mea- 
sured the reflectivity spectrum of ZnCr 2 04 and found a 
large infrared-active (i.r.) phonon splitting (11 cm -1 ) at 
Tn- They invoke a couling between phonons and spins 
to argue that i.r. phonons provide a quantitative measure 
of the spin-Peierls order parameter yj. Yet, a priori, it 
is unclear whether this spin-phonon mechanism can ac- 
count for the size of the effect. 

We study the influence of magnetic order on the 
phonon frequencies of ZnCr204 from first-principles. We 
give evidence that anisotropy induced by AFM ordering 
of spins can account for the large anisotropy measured in 
the i.r. phonons 0. We decouple the macroscopic elastic 
degrees of freedom from those of the spins by calculat- 
ing the phonons of crystallographically cubic ZnCr 2 04. 
This approach is similar to that taken by Massidda et 



al. |fj for MnO and naturally decouples the tetragonality 
induced by the lattice strain at Tjv from that generated 
by the spin pattern alone. Next, in order to understand 
the first-principles results, we revisit the long-standing 
problem of how spins couple to optical phonons Q : 

LO W UJq + A(Sj • Sj) 

We develop an approach where for the first time experi- 
mentally accessible quantities of this theory can be calcu- 
lated from first-principles. Finally, we apply this method 
to calculate the coupling parameter, A, for ZnCr204. We 
find it to be in excellent agreement with that measured, 
providing confirmation of a strong spin-phonon mecha- 
nism responsible for the large phonon splitting. 

First-principles density-functional calculations using 
projector augmented-wave potentials were performed 
within LSDA+U (local spin-density approximation plus 
Hubbard U) |8( as implemented in the Vienna ah initio 
Simulation Package 0, llflj l with a plane wave cutoff of 
500 eV and a 6x6x6 T-centered fc-point mesh. All calcu- 
lations were performed with collinear spins and without 
LS-coupling (inclusion of LS-coupling does not change 
the results). We performed full optimization of the lat- 
tice parameter, a= 8.26 A (exp: a=8.31 A), and anion 
parameter, u = 0.386 (exp: u=0.387), in space group 
Fd3m, where we find excellent agreement with experi- 
ment. The Cr on-site Coulomb, U = 3 eV, and exchange, 
j=0.9 eV, parameters are justified a posteriori from the 
calculation of the magnetic exchange constant as we now 
discuss [llj. 

For zinc chromite, the nearest neighbor (n.n.) exchange 
constant, J, is an order of magnitude greater than all 
next nearest neighbor interactions. J is determined by 
a balance between AFM direct Cr-Cr exchange and FM 
90° Cr-O-Cr superexchange and to a large extent, which 
one wins can be attributed to volume and electronega- 
tivity of the anion. Exchange constants can be ex- 
tracted by mapping first-principles calculations of the to- 
tal energy for different spin configurations at T=0 onto 



a classical Heisenberg model; E s 



(nn) 



Si • S ; 



(sum over z=6 mi's). This gives for the energy of the 
primitive unit cell (4 Cr-ions): Efm=Eq — 24JS 2 and 
EAFM=En + 8JS 2 . From this we calculate J = -2.1 meV 



which compares well with 3 exp = -2.25 meV estimated 
from the measured &cw- 

First- Principles Phonons. —For cubic ZnCr 2 04, group- 
theoretical analysis predicts that the 39 zone-center op- 
tic modes transform according to the following irre- 
ducible representations (irreps) of Oh- IY r . = 4Ti„, 

aman = 3T 2g ®lE g ®lA lg , and T saent = 2A 2u @2E u ® 
2T 2u © ITig, and the three acoustic modes according to 
r acoustic = lTi u . With AFM spins, the ordering pattern 
(the "collinear model" of Ref. 1) induces a tetragonal axis 
chosen here to be z. Under this tetragonal "distortion", 
the point group is lowered from Oh — » D^, and the ir- 
reps previously discussed become reducible, although the 
Raman and i.r. -active irreps still do not mix. Therefore 
we can continue to talk about these modes independently 
and although the Raman modes do split, this splitting is 
small and not our main concern. Still, to demonstrate the 
quality of the phonon calculations for zinc chromite, we 
present our calculated Raman-active frequencies (cm -1 ): 
T 2g : 185 (186), 521 (515), 608 (610); E g : 466 (457); 
Ai g 687 (692), which are in remarkable a gree ment (~1%) 
with experiment (shown in parenthesis) |13| . 

We now proceed to the main part of this letter, the in- 
fluence of magnetic order on the i.r. phonons. Under D^h 
the i.r. modes split according to T\ u — >A 2u ®E u , where the 
1-d A 2u and 2-d E u irreps transform like vectors along 
and perpendicular to the tetragonal axis respectively. 
Additionally, certain rows of the silent T 2u irrep become 
i.r. -active and can in principle mix with those originating 
from T\ u . This mixing is very small and for our purposes 
acceptable to ignore. For the remainder of this letter 
we only discuss the 4 modes originating from T\ u . The 
rows of T\ u do not mix under D^h so that we can still 
talk about the x, y, (E u ) and z {A 2u ) blocks of the dy- 
namical matrix. We compute the i.r. phonons by con- 
structing the three 5x5 blocks (4 i.r. -active, 1 acoustic) 
of the dynamical matrix from Hcllmann-Feynman forces 
on the symmetry-adapted partner functions, f n ,a, where 
n=l,5 and a— x,y,z label the partner function and the 
row of Ti„, respectively. These partner functions trans- 
form like vectors along a under Oh- The partner func- 
tions, fi.a, f 2 ,a, and fA,a, involve atomic displacements 
along d of the entire zinc, chromium, and oxygen sublat- 
tices, respectively. Nontrivial displacements within the 
chromium and oxygen sublattices, whose partner func- 
tions we label fa and respectively, are also possible. 
Although fs j0l and /s ja both transform like vectors along 
d, atomic displacements associated with these partner 
functions take place in the plane perpendicular to the a- 
axis. In Fig. ^ we show fs, z and f^, x which are seen to 
strongly modulate the Cr-Cr bond length. Their impor- 
tance to the present problem will soon become apparent. 

Our calculated i.r. phonons for crystallographically cu- 
bic ZnCr 2 04 with AFM are shown in Table [H labeled ab 
initio. We see that AFM order induces a large tetragonal 
anisotropy in the phonon sector (u} x =to y ^Lu z ) with large 
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FIG. 1: Partner functions f% z (left) and fzx (right) that 
transform like different rows of irrep Ti«. Cr-up spins (white 
spheres) and Cr-down spins (black spheres) are shown while 
oxygen atoms sit at the unoccupied corners are not shown 
for clarity. Arrows indicate direction of atomic displacements 
consistent with f?, a . 

splitting (Alu=lu z -uj x ) occurring for modes (1) and (2), 
Aw 0.12w z , and smaller spitting for (3) and (4), Au = 
0.02cj 2 . (we hold off for now on comparing these results to 
experiment). This mode dependent behavior is reminis- 
cent of the phonon anomalies displayed at T c in the FM 
spinel CdCr2S4 01 where insight was gained by looking 
at the eigendisplacements. We don't repeat the analysis 
hear except to make two points. First, all four modes 
contain Cr-fa partner function (more on this below) so 
that the Cr-Cr bond is highly distorted in the plane per- 
pendicular to a but in (1) and (2) the total distortion 
is such that all Cr-0 bond lengths remain fixed while in 
(3) and (4) the Cr-0 bond lengths, particularly along a, 
are strongly modulated (we will see what this implies). 
Second, the presence or absence of Zn-/i has no effect 
on the spin-phonon coupling in contrast to what is com- 
monly proposed. We understand this from our calcula- 
tions of the q=0 component of the force constant matrix, 
C n ,a;n>,ot>, where we find Ci iX; i jX = Ci. z -i. Zl as you would 
expect for cubic symmetry. We found similar results in 
CdCr 2 S4, CdCr 2 Se4, and HgCr 2 Se4 which suggests that 
the previous interpretation of phonon anomalies in these 
FM spinels should be revisted 0. 

Our first-principles calculations clearly show that AFM 
spin-ordering induces a large i.r. phonon splitting as 
shown in Table [I] Sushkov argued that the large phonon 
splitting of Ti u (2) was caused through modulation of di- 
rect Cr-Cr exchange by fa, and since T lu {2) contained the 
largest amount of this partner function, it was the only 
splitting measurable. In Fig. ^ the anisotropy imposed 
by the magnetic order on this mode is clear. Supporting 
their argument is our calculations of the q—0 component 
of the force constant matrix. Looking at the diagonal 
components we find that the anisotropy, C n ^ n ^ z -C n ^ n%x 
(which equals zero for a cubic system), is «20% for / 3 , 
being 3 times larger than that due to any other partner 
function. Second, Ti„(2) contains twice as much f% com- 



3 



TABLE I: Infrared-active, Ti u , phonons frequencies, cm 
of ZnCr 2 4 . 



Experiment Present Theory 
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ab intio 




13K 


9K 


PM 


AFM 


Jl 




AFM 






x z 




X z 






x z 


(1) 


186 


186+5 


189 


186 191 - 


■240 





174 198 


(2) 


371 


368 379 


366 


361 372 - 


-999 


22 


342 392 


(3) 


501 


501 


514 


514 517 - 


■254 


225 


510 526 


(4) 


619 


619 


621 


621 623 


-66 


257 


620 630 



pared with the other i.r. modes. But in contrast to what 
is seen experimentally, the phonon frequencies of all four 
modes are split. Additionally, the calculated splitting of 
Ti u (2), Aw=50 cm -1 , is about 4.5 times larger than the 
measured splitting, Acj=11 cm . One question then 
becomes how do we interpret our phonon calculations to 
be consistent with experiment? Also what do they teach 
us about the spin-phonon mechanism and the role of di- 
rect exchange? Clearly a more systematic approach to 
understand the first-principles calculations is desirable. 

First- Principles Model— We begin by writing the total 
energy, E=Eo+E p h+E sp in, of a system of i.r. phonons, 
^'p' i= 2 r l = { n i OL } and C rj ^ rj '=8 aa iC n ^ l i , 

and spins, E sp i„=— JijSi-Sj. Without the spin-orbit 
interaction the spin and phonon sectors are still indirectly 
coupled due to a positional dependence of the exchange 
interactions 0. This dependence is rather complicated 
because in addition to the positions of the two spins, the 
positions of the surrounding anions are important for the 
superexchange process, and may even influence direct ex- 
change In principle Jij depends on the positions of 
all N magnetic and non-magnetic ions, Jij(r\,r2, ■ ■■?"./v)- 
A solution then becomes clearer if we write J in terms 
of f n since these functions provide a complete basis to 
represent all possible i.r. atomic motion. For small dis- 
tortions, J can be expanded in \f n \ leading to 

d 2 J ■ 

to lowest order in SfSj. In general, symmetry of the 
magnetic system can reduce this expression further. For 
spinels, D±h does not induce couplings between differ- 
ent rows of O h , i.e. C VtV >=C na ^ a = C n ^(a), so in 
fact we only need to calculate derivatives of the type, 
d 2 J/df n df n > (magnetic order other than D 4 h could in- 
troduce derivatives connecting different rows of Oh but 
will not change the value of the intra-row type). Now 
since the crystallographic group is cubic there are only 
two distinct changes in Jy depending on whether the 
spins lie in a plane perpendicular to a or not. With this 
let us define d 2 J l j/df v df v > = J" n n , ^ V fy • a = 



and d 2 J lj /df ri df ri > = J'J, ^ V hj • a ^ where 
is the direction vector linking nearest neighbor spins, Si 
and Sj. Then we can write 

C n>n '(a) = C n ^(a)-J'l n n , ^(S r Sj)-J||' n>n , ^(S r Sj) 

fx ni 

with the first sum over the two neighbors in a plane per- 
pendicular to a and the second over the remaining four. 

We have considerably reduced the problem by writ- 
ing the exchange constants in terms of the partner func- 
tions and by applying symmetry, but we are still faced 
with the challenge to calculate J" . For this we propose a 
procedure that is similar to the approach of calculating 
exchange constants - by mapping onto a Heisenberg-like 
model. As such let us write 







— 8 J'Ln ,ra' 


.s 2 






^PM 


OJ _Ln,n' 


■ S 2 


+ 16J||„,„^ 2 




n PM 
— °n,n' 


+ & J ln,n> 


■S 2 





where C and (J are the force constant matrices 
calculated from first-principles with FM and AFM or- 
der respectively. Here we have three equations and three 
unknowns. In Table [I] we show the T—0 PM phonon 
frequencies extracted from the model which are found 
to be in excellent agreement (1-2%) with those mea- 
sured at 13K. To test the validity of the model w.r.t. 
first-principles we perform additional phonon calcula- 
tions with different magnetic order, for example, ferri- 
magnetic (FiM) order with one Cr-ion spin-down, the 
other three spin-up. For this spin configuration accord- 
ing to our model C FiM {x)=C FiM (y) =C FiM (z)=C PM . 
If we now perform first-principles calculations of the 
phonons with FiM magnetic order, we find that the fre- 
quencies are in exact agreement (<lcm~ 1 ) with those we 
extracted from the model, i.e. uj FiM — lu pm 15]. As a 
further test we considered additional FiM spin configu- 
rations within the conventional unit cell (eight formula 
units). The model reproduced the first-principles calcu- 
lations to within 1 cm -1 in every case. 

Discussion.— Having established the model correctly 
captures the influence of magnetic order on the first- 
principles phonons, we turn our attention to calculat- 
ing experimentally accessible quantities of the problem. 
First, on the role of direct exchange, it is telling that 
we find d 2 J ± /df 3 df 3 w-0.172 meV/A 2 , more than -10 
times larger than any other component of </" or J'J. This 
anomalously large value is consistent with the anisotropy 
in the force constant matrix and, as we now show, as- 
sumed correctly by Sushkov to originate from the modu- 
lation of direct exchange, J^. For small displacements 
Jd goes like Jd{R c )~Jd^ a ^ Rc , where R c is the dis- 
tance between neighboring Cr-atoms. Assuming the 
contribution from superexchange (SE) is small, we find 
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a=9.05A _1 which compares well with the measured value 
of o^S.OA^ 1 . This must be considered somewhat fortu- 
itous since it is not clear that either the theoretical or the 
experimental approach isolates direct exchange. What 
does suggest the dominate role of direct exchange in J" 3 3 
are our calculations for CdCr 2 04. This compound has a 
larger lattice constant and less negative J than ZnCr2 04 
so that direct exchange is expected to be weaker. Repeat- 
ing the phonon calculations and mapping procedure for 
CdCr 2 4 , we find d 2 J^/df 3 df 3 « -0.09 meV/A 2 . Con- 
sidering a small contribution from SE is expected this 
compares well with a 2 J=-0.07 meV/A 2 . Furthermore, 
no other components of J" or J 'J change appreciatively 
from those of ZnCr 2 04, consistent with the fact that the 
only partner function to change R c is f 3 . 

The relation between J" 3 3 and direct exchange being 
established, we return to the discussion of the phonon 
splitting in ZnCr 2 Oi. To quantify the mode dependence 
and compare with experiment, let J x —u\j"u\, where 
ma are the PM dynamical matrix real-space eigenvectors. 
Then we can write for the phonon anisotropy 

uj x (z) - uj x (x) « 4 — Si • S a - Si ■ S 4 ) 

where (Si • S 2 — Si • S4} is defined as the spin-Peierls order 
parameter, f sp , In Tabled we see that the anisotropy 
experienced by T\ u (1) and (2) is due to J", which is 
dominated by J" 3 3 and subsequently direct exchange. In 
contrast, J7j" plays a large role in the anisotropy of T\ u (3) 
and (4) . This is due to the strong modulation of the Cr- 
O bonds which we suspect modifies FM superexchange 
and may also explain the sign. We can now compute 
the spin-phonon coupling constant, A, for T\ u (2). We 
find A 2 =4(J|' — J")/ljpm~11 cm -1 which agrees well 
with experiment, A 2 = 6-10 cm -1 . So from calculations 
of just the phonons for selected spin configurations, we 
are able to extract the spin-phonon coupling parameter, 
X, from first-principles. It is also now clear that the cal- 
culated splittings are significantly larger than experiment 
because the spin-Peierls OP for collinear AFM, f sp = 4.5, 
is significantly larger than experiments (extracted values 
range from f sp = 1.1 to 1.8). In Fig. we show how 
the splitting for T\ u (2) evolves with the spin-Peierls OP, 
where at f sp = 1 we recover the experimentally observed 
splitting. If we substitute this value of the spin-Peierls 
OP into our model, we find that the splittings of the other 
modes are significantly reduced (2-5 cm -1 ) as shown in 
Table H] Although additional effects of the structural 
transition at T^r still need to be worked out, perhaps 
such small splittings are not easily measurable. But, it is 
clear, if ZnCr 2 04 could be prepared in a state where f sp 
attained its maximum value the splitting would be 2.5 to 
4 times greater than that currently measured ||. 

In summary, we have developed a method to calculate 
parameters of a spin-phonon theory from first-principles. 




FIG. 2: Phonon frequency of Wa(-S) and U2(x) as a function 
of the Spin-Peierls order parameter (Si • S2 — Si • S4). 

The method accounts for direct exchange which was 
shown to be responsible for the large magnetic induced 
anisotropy in the phonon channel of AFM ZnCr 2 04 but 
also equally for superexchange. As a consequence a nat- 
ural solution to the problem of phonon anomalies in FM 
spinels is also revealed as a balance between these two 
processes. Q We anticipate that the approach developed 
will also find application to magnetoelectric problems 
where spin fluctuations coupled to optical phonons have 
been proposed to explain the magnetocapacitance effect. 
Given the breath of such problems today, how the spins 
couple to the lattice and an approach to calculate the 
relevant parameters of such a theory from first-principles 
are important questions that needed to be addressed. 
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